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Muldmode spinless p-wave superconducting wires with a width W much smaller than the superconducting 
coherence length ^ are known to have multiple low-energy subgap states localized near the wire's ends. Here 
we compare the typical energies of such endstates for various terminations of the wire: A superconducting wire 
coupled to a normal-metal stub, a weakly disordered superconductor wire, and a wire with smooth confinement. 
Depending on the termination, we find that the energies of the subgap states can be higher or lower than for the 
case of a rectangular wire with hard-wall boundaries. 
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I. INTRODUCTION 

In the current search for Majorana fermions in nano-wire 
geometries' ^ an important theoretical challenge is to under- 
stand the multiplicity of possible fermionic bound states that 
can form at the ends of the wire and how a possible Majorana 
bound state can be identified among them. This is particu- 
larly relevant for multichannel geometries, in which fermionic 
states localized near the ends of the wire are expected to oc- 
cur at energies much smaller than the excitation gap for bulk 
excitations if the wire width is much smaller than the super- 
conducting coherence length. In this article we explore the 
dependence of these sub-gap end-states on the details of the 
termination of the wire and on impurity scattering. 

The interest in isolating Majorana fermions arises because 
their non-local properties and non-abelian braiding statis- 
tics render them potentially useful for fault tolerant quantum 
computation. '' '' Majorana fermions occur — at least theo- 
retically — at the ends of one-dimensional spinless p-wave 
superconductors.'*^ Recent proposals suggest ways of engi- 
neering solid-state systems that effectively behave as spinless 
p-wave superconductors by combi ning an s-wave supercon- 
ductor and a topological insulator," '^ a semiconductoJSE^ 
or ferromagnet.'^ ^' Building on the proposals of Refs.fTSjand 
[T6l two experimental groups have reported an enhanced tun- 
neling density of states at the ends of InAs and InSb wires 
in proximity to a superconductor, consistent with the exis- 
tence of Majorana bound states at the ends of these wireSjI^^SI 
whereas a number of other groups claim the observation of 
Majorana bound states using different methods.l^^l^Il 

Whereas the original proposals for Majorana fermions in 
wire geometries focused on one-dimensional systems, it is 
by now well established that the topological superconduct- 
ing phase with Majorana end states may persist in a quasi- 
one-dimensional multichannel setting. ^^ "'^ A difference be- 
tween the quasi-one-dimensional and one-dimensional set- 
tings is, however, that a possible zero-energy Majorana state 
localized at the wire's end may co-exist with other fermionic 
sub-gap states, analogous to those found in vortex cores of 
bulk superconductors.^^** For the case of an iV-channel spin- 
less p + ip superconductor with a rectangular geometry and 
with width W much smaller than the superconducting coher- 
ence length ^, three of us recently showed that the number of 



such fermionic subgap states is ^ N/2, and that their typical 
energy is ^typ ~ ^ being the superconducting gap 

size.^^ The lowest-lying and highest-lying fermionic subgap 
states have energies Smin ^ £typ/N In N and Smax ^ Netyp, 
respectively. The fermionic subgap states also exist in a non- 
topological phase without zero-energy Majorana end-state, 
thus posing a potential obstacle for the identification of the 
topological phase through the observation of an enhanced 
density of states near zero energy. 

In a recent article. Potter and observe that the de- 

pendence of the energy of the lowest-lying fermionic subgap 
state on system parameters changes qualitatively if the rect- 
angular geometry of Ref. 21 is replaced by a geometry with 
rounded ends. They point out that the calculation of the en- 
ergy of the fermionic subgap state for the rectangular geome- 
try is plagued by a subtle cancellation, which does not appear 
for a generic wire ending. In particular it was found in Ref.l28l 
that the lowest-lying fermionic subgap state has an energy sig- 
nificantly above the prediction of Ref. 27 for a wire with width 
W ^ £, and rounded ends. 

Motivated by these observations we present here a detailed 
investigation of the effect that the wire termination has on the 
energies of the fermionic subgap states for the multichannel 
spinless p + ip superconductor Remarkably, we find that, 
depending on the details of the wire ending, the energies of 
the fermionic subgap states can be significantly above, as well 
as below the rectangular- wire case of Ref. 27 We find an 
increase of the energies of the subgap states if an arbitrarily- 
shaped normal layer is attached to the wire's end, the magni- 
tude of the increase being consistent with the estimate of Ref. 
l28lfor a wire with rounded ends. On the other hand, the pres- 
ence of im purities — weak enough to preserve the topological 
phasPE- - on average reduces the energies of the fermionic 
end states below the estimate of Ref. 27 , while a smooth con- 
finement (with a slowly increasing potential energy providing 
the confinement along the wire's axis) leads to even smaller 
energies of the fermionic subgap states. 

Our results are derived for the two-dimensional spinless 
p + ip superconducting strip of width W. The model of a 
spinless p + ip superconductor is an effective low-energy de- 
scription for the various proposals to realize one-dimensional 
or quasi-one-dimensional topological superconductors, pro- 
vided the number N of propagating channels at the Fermi 
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level is chosen equal to the number of spinless (i.e., helical 
or spin-polarized) channels in the case of the semiconductor 
or ferromagnet proposals. (The edges of a topological insula- 
tor always have = 1, so that a multichannel p + ip model 
is not relevant in that case.) A mapping between the spinless 
p-wave model and the semiconductor-wire proposals is given 
in the appendix. 

The remainder of this article is organized as follows: In 
Sec. [n] we briefly review the symmetries of the model (j2|i and 
the reason for the appearance of multiple low-lying states if 
the wire width W is much smaller than the superconducting 



coherence length ^. In Sec. Ill we describe a scattering theory 
of fermionic subgap states with arbitrary wire endings. Sec- 
tion |IV] discusses the p + ip model with weak disorder, while 
the effect of a smooth potential at the wire's end is discussed 
in Sec.lv] We conclude in Sec.lVll 



II. p + ip MODEL 

Our calculations are performed for a two-dimensional spin- 
less p + ip superconductor, which is described by the two- 
component Bogoliubov-de Gennes Hamiltonian, which we 
write as 



H — Hq + Hy + Hy, 



(1) 



with 



2m 



Hy = -A'pyTy, 



(2) 



Here t^, Ty, and are Pauli matrices in particle-hole space, 
A' specifies the p-wave superconducting order parameter, 
/Lt = h'^kp/2m and m are the chemical potential and elec- 
tron mass, and V{r) a potential that describes the confine- 
ment at the ends of the wire as well as the scattering off im- 
purities. The two-dimensional coordinate r — {x, y), where 
< y < W, with hard-wall boundary conditions at y = 
and y = W. The superconducting order parameter derives 
from proximity coupling to a bulk superconductor, so that no 
self-consistency condition for A' needs to be employed. 

Hypothetical end states are localized within a distance 
of the order of the superconducting coherence length ^ = 
h{A'm)~^ from the wire's ends. For thin wires with W <^ 
it is a good starting point to analyze the Hamiltonian H = 
Hq + Hy without the term Hy. The Hamiltonian Hq has a 
chiral symmetiy^ TyHoTy = —Hq, and there exist 
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FIG. 1. (Color online) Schematic picture of the spectrum of low- 
energy excitations of the a. p + ip wire as a function of its width 
W. The gap for bulk excitations closes at those values of W for 
which (W/n)^kp — is an integer. When the bulk gap is fi- 
nite, there are low-energy subgap states localized near the ends of 
the wire. In the text, we use Eniin to denote the energy of the lowest- 
lying fermionic subgap state, £typ for the typical energy of a subgap 
state, and Emax for the energy of the highest-lying fermionic subgap 
state. 



Majorana end states, since Hy preserves the chiral symmetry, 
although it may change the boundaries of the phases with dif- 
ferent N if Hy is nonzero in the bulk of the wire. In contrast, 
the term Hy breaks the chiral symmetry and couples the N 
Majorana bound states, giving rise to (generically) int {N/2) 
fermionic states at each end and a single Majorana end state 
if N is odd. If ^ ^ the splitting of the end states is small 
in comparison to the bulk energy gap A = A'hkp, and the 
resulting fermionic states cluster near zero energy. 

A schematic picture of the end-state spectrum as a func- 
tion of W is shown in Fig. [T] The end states are character- 
ized by the energy Enun of the lowest-lying fermionic end 
state, the typical end-state energy £typ, and the energy Emax 
of the highest-lying end state. For small N these three energy 
scales are comparable, but for large they may differ consid- 
erably. The energy Emin serves as the "energy gap" protecting 
the topological state and sets the required energy resolution if 
the presence or absence of a Majorana end state is detected 
through a tunneling density of states measurement. 

The specific case of a rectangular wire geometry, with hard- 
wall boundary conditions at each end of the wire and without 
disorder, was investigated in Ref.lZT] We now investigate two 
other possible terminations, as well as the effect of disorder 
on the energies of subgap endstates in multichannel spinless 
p-wave superconducting wires. 



N = mt[{W/TT)^ 



(3) 



Majorana bound states at each end of the wire.lSiniThe step- 
wise increase of the number of Majorana end states for wire 
widths W such that (W/n) \/kp — is an integer is accom- 
panied by a closing of the bulk excitation gap of Hq. Inclusion 
of the potential term Hy does not lift the degeneracy of the 



III. NORMAL-METAL STUB 

In this section, we consider a quasi-one-dimensional spin- 
less p + ip superconductor without disorder and coupled to 
a normal-metal stub at its end. We choose coordinates, such 
that the spinless superconductor occupies the space x > 0, 
< y < W, see Fig. |2] Such a wire ending is relevant, e.g., 
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for the experimental geometry of Ref.|22l in which a topolog- 
ical phase is induced in a semiconductor nanowire by laterally 
coupling it to a superconductor, while a part of the wire sticks 
out from under the superconductor and is pinched off by a gate 
at a finite distance. 

We take the Hamiltonian of the normal stub to be real and 
symmetric, in order to preserve the chiral symmetry of the 
Hamiltonian Hq. Following Ref. 27 we first solve for the 
wavefunctions i/'*^-'^ of the N Majorana modes for the Hamil- 
tonian Hq and then treat Hy in perturbation theory. The po- 
tential term Hy is set to zero throughout this calculation. 

The Majorana states have support in the normal stub as well 
as in a segment of the superconducting wire of length ^ ^. In 
the superconducting region a; > the wavefunctions tp of the 
Majorana states can be written as 



i-(f>n-ir) + a„+(/)„+(r), 



(4) 



where the basis states (f>n±, n — 1,2, . . . , N, read 

pj7r/4 



4'n±ir) 



with 



-i7r/4 



whk„ \ w J 



- - {nir/Wy 



(5) 



(6) 



The basis states (/)„± have been normalized to unit flux. The 
above expressions for the basis states and their normalization 
are valid up to corrections of order (W/^)^, which we neglect 
throughout this calculation. 

The coupling to the normal-metal stub imposes boundary 
conditions on the coefficients a„±, which we express in terms 
of the scattering matrix Snn' of the normal stub. 



V 9* 



(7) 



Because the Hamiltonian of the normal stub is real and sym- 
metric, the scattering matrix Snn' is unitary and symmetric, 
Snn' — Sn'n, which cnsurcs that the 2N equations (|7]i have 
TV independent solutions, corresponding to the N Majorana 
end states. 

For finding an explicit representation of the N Majorana 
states ^^^\ j = 1,2, . . . , N wc use the fact that the scat- 
tering matrix S and the Wigner-Smith time-delay matrij(E2E2] 
Q — ihS^dS/dfi of the normal stub can be simultaneously 
decomposed as 



S = U^U, Q = C/^diag (n, . . . , tn)U, 



(8) 



where U is an N x N unitary matrix and the Ti > 0, i — 
1,2, . . . , N, are the so-called "proper time delays". With this 
decomposition, a solution to the boundary conditions (j?]) is 
given by 



Un 



U* 



J 



1,2,.. 



(9) 




FIG. 2. (Color online) Schematic drawing of a spinless p-wave su- 
perconducting wire (S) coupled to a normal-metal (N) stub at one 
end. The top panel shows a rectangular stub, the bottom panel shows 
a chaotic cavity attached to the superconducting wire. 



The N states that are defined through these coefficients. 



(/>„_(r) + a|/j0„+(r), 



(10) 



are Majorana modes (they satisfy = t^^^^''), but they 

are not necessarily orthonormal. In order to construct an or- 
thonormal set, we first calculate the scalar product Mji of the 
modes ^ 



(11) 



' Stub 

E2m^ 



Re 



Tl -U* 
TT .TJ* A. ^"J^"' 



- 2t,-(5 



Here we used the relation between the Wigner-Smith time de- 
lay matrix and the normalization of scattering states in order 
to perform the integration over the sub, see Ref. 43 The 
overlap matrix M is real, positive definite, and symmetric. 
It is manifestly diagonal if the scattering matrix S and the 
time-delay matrix Q are diagonal or in the "large-stub limit", 
which is defined as the limit in which the mean inverse dwell 
time h/f is much smaller than the superconducting gap. In 
both cases, one obtains an orthonormal basis for the N Ma- 
jorana modes by setting tp^J^^ = tp'-^'> / ^ Mjj . In the general 
case, M is not diagonal, however, and one has to construct 
an orthonormal system with the help of the orthogonal trans- 
formation O that diagonalizes M, i.e., O^MO — X^, where 
A = diag (Ai, A2, . . . , Xn) is a diagonal matrix with positive 
elements. The corresponding orthonormal basis set one thus 
obtains reads 



N 

{r) = J24^^"\r)OnjK. 



(12) 



Inclusion of Hy, which breaks the chiral symmetry, gives 
rise to a splitting of the N degenerate Majorana end states 
constructed above. With respect to the unnormalized states 
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this splitting is described by the N x N matrix 



H' 



(1) 



4iA'TO ^ nn'[l - (-1)"+"'] 



W 



7in 



E 

± I- 



Im 



^ {n'^ — in? ) \/ kn h 



,(13) 



where we neglected corrections smaller by a factor of order 
(W/^)^. The matrix Hj^'' is antisymmetric and purely imag- 
inary, which ensures the existence of a single zero-energy 
bound state if N is odd. In order to find a true effective Hamil- 
tonian H^^\ the eigenvalues of which represent the energies 
of the fermionic end states, one should transform to the basis 
of orthogonal states t/;'^^ introduced in Eq. (12\ , 



A A 



(14) 



In the special case N — 2, this transformation can be carried 
out for an arbitrary scattering matrix S and the energy of the 
resulting single fermionic bound state is 



\H- 



(1)1 

12 I 



22 



(15) 



We now discuss two particular realizations of a metal stub in 
detail. 




FIG. 3. (Color online) Typical and maximal energies of fermionic 
subgap states in a spinless p-wave superconductor with a rectangular 
normal-metal stub of length L as a function of L/W for different 
channels numbers (A'^ = 15, 27, 55, 99). The maximal energies 
Eniax have a finite- A*' correction of order Sm/V^, which is why the 
curves for Emax still show an N dependence for large A''. 




A. Rectangular stub 

First, we consider a rectangular stub of length L attached 
to the spinless p-wave superconducting wire, see Fig.|2^. For 
this geometry, both the scattering matrix S and the Wigner- 
Smith time-delay matrix Q are diagonal. 



J2ik„L s 



n - 2mL 

flKri 



(16) 



with kn given by Eq. (|6]l. Since there is no mixing between 
different channels, the zero energy modes V^'^^ already form 
an orthogonal basis. The effective Hamiltonian in the normal- 

= if j + Hs even 

sm[L{kj ± ki)] 



(1) 



ized basis V'^' = V'^^'^ / has Hji 



and 



Wi^ + 2L){P - P 
2W cos[L{kj ± ki 

^ ^ {ki± kj Y 



E 



ki ± ki 



(17) 



if j + I is odd, up to corrections smaller by a factor of order 

The second term in the effective Hamiltonian ( [T7] i is smaller 
than the first one by a factor of order W/ ^. However, only this 
second term contributes in the limit i = in which there is 
no normal metal stub.^^ This is a variation of the cancella- 
tion effect pointed out by Potter and Lee.^** We now analyze 



FIG. 4. (Color online) Level density of fermionic subgap states for 
a rectangular stub in the limit of large A*', for L/W = 0.1 (left) 
and L/W — 3 (right). The level density is measured in units of 



the eigenvalues of the effective Hamiltonian H^^'' for finite L, 
when the first term between brackets dominates. 

Since no closed-form expressions for the eigenvalues of 
H'^^'> could be obtained, we numerically diagonalized H^-^^^ 
and investigated the dependence of the minimum, typical, and 
maximal positive eigenvalues on the ratio L/W as well as the 
channel number N. For L ^ W, this analysis gives 



yp 



WA 



(18) 



see Fig. [3] The maximum and minimum energies of the sub- 
gap states scale as Emin ~ £typ/N, e^ax ^ £typ- A similar 
analysis for fcp ^ ^ i <C gives estimates for Emin, etyp, 
and £,nax which are smaller by a factor L/W, whereas for 
L <C fcp ^, they are smaller by a factor L^kp /W. A crossover 
to the results of Ref. f^Tl takes place for L < (VF^^fO^^^- 
In the limit of large N the energies of the fermionic subgap 
states are best described through their level density, which is 
shown in Fig.|4] 
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FIG. 5. (Color online) Level density of fermionic subgap states in 
the small-cavity limit (left) and large-cavity limit (right). The level 
density is measured in units of — N/es and z^i = N/ei for the left 
and right panels, respectively. 



B. Chaotic Cavity 

As a second example, we consider a chaotic cavity attached 
to the end of the superconducting wire, see Fig. |2]5. In this 
case, the unitary matrix U is randomly distributed in the uni- 
tary group,"^ whereas the proper delay times have the proba- 
bility distributiorPSl 



N 



P(ri,...,Tjv)(x[]^(^j)' 



_-3N/2-l Nf/2Tj 



i<j 



(19) 



with the average delay time f. In this case, the matrix U 
is not diagonal, and the prescription of Eq. ^\A\ needs to be 
used in order to construct the effective Hamiltonian H^^"^ for 
the low-energy subgap states. As in the previous case, we 
could not obtain closed-form expressions for the energies of 
the fermionic subgap states and had to resort to a numerical 
analysis, in which the unitary matrices U were generated ac- 
cording to the Haar measure on the unitary group and the time- 
delays Ti according to the probability distribution given above, 
following the method described in Ref. 46 This analysis gives 
different results for the limiting cases of a "small cavity" and 
a "large cavity", corresponding to the inverse mean dwell time 
h/f large or small in comparison to the superconducting gap 
A. 

Small-cavity limit. In the small-cavity limit, the normaliza- 
tion of the N Majorana states t/;^-'-' is dominated by the inte- 
gration over the superconducting wire. Not counting the Ma- 
jorana states, the excitation spectrum of the cavity has a gap 
comparable to the bulk excitation gap A. Upon including Hy 
one obtains N fermionic subgap states, which have a typical 
energy 



yp 



(20) 



and Eniax ~ Etyp, Emin ~ £typ/N. The prccisc location of the 
subgap states depends on the precise scattering matrix of the 
cavity. The mean level density for an ensemble of cavities is 
shown in the left panel Fig. |5] 



Large-cavity limit. In the large-cavity limit, the overlap ma- 
trix Mji is dominated by the in-cavity parts of the wavefunc- 
tions, so that the Majorana modes ■i/;^''^ are already orthogonal 
and the effective Hamiltonian Hjj^ = Hjj^/ y^ArjTi, with 



Eq. 



Not counting the Majorana states, the 



H^i given in 

cavity's excitation spectrum has a gap of order = h/nf^ 
where f is the mean dwell time in the cavity. In this case, the 
typical energy of the fermionic subgap states is 



Etyp ~ ei 



(21) 



while Emax ^ Etyp and £niin ^ Etyp/N. The mean level den- 
sity of the subgap states for an ensemble of cavities is shown 
in the right panel of Fig. [5] 

IV. p + ip MODEL WITH DISORDER 

Whereas strong disorder is known to destroy the topologi- 
cal superconducting phase in the p + ip model in one dimen- 
sion, weak disorder with mean free path Z > preserves 
the topological phase. In this section we investigate the ef- 
fect of weak disorder on the energies of the fermionic subgap 
states in a multichannel rectangular p + ip model. Because 
the disorder is necessarily weak (strong disorder suppresses 
the topological phase), the effect of disorder can be treated in 
perturbation theory. 

Starting point of our analysis is the chiral-symmetric 
Hamiltonian Hq, which has N normalized Majorana bound 
states I t/i'^'), j — 1,2, ... ,N at each end of the wire. We 
take a rectangular geometry, with a wire end and hard-wall 
boundary conditions at a; = 0, and take the potential V{r) to 
be a Gaussian white noise potential with mean {V{r)) = 
and variance 



{Vir)Vir')) = ^^Sir~r'), 



(22) 



where / is the mean free path and vp = hkp/m the Fermi 
velocity. In our perturbative analysis we treat both the impu- 
rity potential V and the transverse superconducting order as 
perturbations and write 



H = Ho + U, 



(23) 



where U = Hy + Hy contains the superconducting correla- 
tions coupling to py as well as the impurity potential. 

The effective Hamiltonian Hcs describing the splitting of 
the N Majorana states into fermionic subgap states can be 
found using the degenerate perturbation theory of Kato"*^ and 
Bloch."*^ (For additional details on this methodology see also 
Refs.EniandEII) Defining P = E I I as the pro- 

jector onto the zero-energy subspace and Q = 1 — P, we can 
then write using Bloch's method 

Q Q Q 

H^s = PUP - PU-^UP + PU-^U-^UP 

Ho Hq Hq 



- ( PU-%UPUP + PUPU-%UP 

2 V iJg ^^n 



■(24) 
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It is essential to note that the disorder potential V{r) alone 
cannot lift the degeneracy of the Majorana end states at any 
order of the perturbation theory. This can be understood di- 
rectly from the observation that the disorder potential V{r) 
does not break the chiral symmetry of the unperturbed Hamil- 
tonian Hq that is responsible for the A^-fold degeneracy. On 
the level of perturbation theory this can be understood im- 
mediately through the particle-hole symmetry present in the 
Majorana bound states and the knowledge that for each per- 
turbative diagram that connects Majoranas through the posi- 
tive energy bulk states there is a cancelling path through the 
negative energy states. 

Keeping terms to first order in Hy and up to second order 
in Hy only, we obtain 



=i7(l)+i/(2)+i/(3a)_^(3b)^ 



(25) 



with 

Q 



H 



(2) 



= (V-W) J- ff^. J- + permutations, 



(3b) 



where 



k 



(2) 



3 1 



^0 



(0\ 



(26) 



(27) 



The effective Hamiltonian H^s is antisymmetric, which im- 
plies that the diagonal elements of all the above terms are zero. 
The first-order term H^^^ describes how the transverse super- 
conducting correlations lift the degeneracy of the N Majorana 
modes in the absence of disorder. The second-order term iJ 
is linear in the disorder potential. Its elements are random 
variables with zero mean and standard deviation that does not 
appreciably change with ^. The third order term contains two 
terms, the first of which is also a random variable with zero 
mean and with a root-mean-square proportional 

The term H^^^^ contains corrections to the effec- 
tive Hamiltonian arising from the renormalization and re- 
orthogonalization of wavefunctions at the first order of the 
perturbation theory. Since this term is a weighted sum of first 
order elements H^^\ it is the only one of the higher-order 
perturbation corrections that gives a systematic dependence of 
energies on the disorder strengths. To see this in more detail, 
it is instructive to separate the diagonal and the off-diagonal 
elements of V^^'^ in the expression 

fori/(3b) 



H 



(3b) _ 



(2)ff(l) 



+ \T^i^H^+H^yl^). (28) 

The first term here is the most important because the weights 
Vj^i^ are positive definite random variables. A simple scal- 
ing analysis predicts that these variables have both mean and 



standard deviation proportional to the ratio S^/l of coherence 
length and mean free path. This term effectively renormalizes 
the entire first order contribution, on average driving the ener- 
gies of the fermionic subgap states towards zero. The second 
term, which contains the contribution from the off-diagonal 
elements of V'^'^\ is less important because the disorder po- 
tential here connects different Majorana modes. These matrix 
elements are therefore randomly distributed with zero mean 
and a root mean square proportional to the coherence length. 

Motivated by these observations, we write the effective 
Hamiltonian in the form 



(29) 



where c = (l/N^) V^^^ is a number of order unity, and 

H' = ^(2) + (^(3a) _ l|y(2) _ (30) 

The correction H' has zero mean. 

We have numerically diagonalized a lattice version of the 
Hamiltonian Q in order to provide numerical evidence for 
the applicability of the above arguments. Details of the re- 
lationship between the continuum and lattice models can be 



I Results of the numerical simulations are 
For weak disorder, the perturbation ij'^^ 



found in Ref. 
shown in Fig. 

dominates the response of the fermionic subgap states, and 
the energies of the fermionic subgap states may both increase 
or decrease, depending on the specific realization of the dis- 
order potential. While large fluctuations persist, for stronger 
disorder the quadratic-in-disorder perturbation (^b) jg^^jg iq 
a systematic decrease of the energies of the fermionic subgap 
states, which is well described by a linear dependence on ^//, 



consistent with the first term in Eq. ( 29 1 



V. SMOOTH POTENTIAL AT WIRE'S END 

In this section we consider a wire which is terminated by 
a smooth potential V{x), as shown schematically in the in- 
set of Fig. [T] In order to address this scenario we solve the 
Bogoliubov-de Gennes Hamiltonian in the WKB approxima- 
tion. Without the transverse pairing term Hy there are N Ma- 
jorana states t/;^-'' with wavefunction 



g»7r/4 
p — i7r/4 



Xj(a;)sin('^) , (31) 



where the functions Xj{x) take the form 

Xjix) 



e^^'? cos[$j(a;)] 



if a; < x-i 



if a; > x-j 



(32) 
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FIG. 6. (Color online) Distribution of energies of fermionic subgap 
end states in a spinless p-wave superconductor with N = 7 chan- 
nels (dots), as a function of ^/L For small amounts of disorder the 
term ff'^' dominates, pushing the subgap energies up or down with 
equal probability. At stronger disorder the term in oc ^/l 

eventually dominates and pulls all energy levels towards zero. The 
red lines indicate the mean calculated from the local distribution of 
eigenvalues. The black lines, which are a linear fit to the mean val- 
ues in red, share an approximate common intercept at the horizontal 
axis at 5// = — 2.2. Dotted grey lines indicate the unperturbed 
energies. Energies are measured in units of e[yp = AM^^/^^. The 
lattice parameters used in the numerical calculation correspond to 
kpW ^ 23andkF^ ~ 320. 



where $j = n/A — dx'kj{x'). 



kj (x) = J 2m fj, — V{x) — 



2mW^ J C 



n,{x)^J--2m ^^,-V{x)-^^ 



ilj is the normahzation constant, and Xj is the transverse- 
mode-dependent classical turning point, defined as the solu- 
tion of kj{xj) — 0. Inclusion of the transverse pairing Hamil- 
tonian Hy lifts the degeneracy of the N zero energy Majorana 
end states, where the energy splitting is given by the eigen- 
values of the antisymmetric matrix with elements H^^'' = if 
j — I even and 



(i)_16A' Ij (1) 



X 



(2) 



X 



(3)^ 



(33) 



if i — I odd, with, for j < I, 

— 2x/^~f^ dx' Kj(x')~ J^^ dx' Ki{x') 

dx 













J — oo 






xf = 




J Xj 






xf . 






J Xi 



dx 



dx 



■'^i cos[$j(a:)] 
2^Vlj^ikjK,i 
g-22;/c cos[<i>j(x)] cos[$i(a::)] 



^/rijflikjki 



(34) 



FIG. 7. (Color online) Energies of fermionic subgap end states, 
in a 6-channel p + ip wire as function of the adiabaticity parame- 
ter a (red). The flat black lines indicate the energies for hard-wall 
boundary conditions (a — 0). Energies are measured in units of 
e[yp — AW^ /^'^. The lattice parameters used in the numerical cal- 
culation correspond to kpW ~ 19 and kp£^ ~ 90, a = Bfi. The 
inset shows the potential profile used in the calculations. 



Figure |7] shows numerical calculations for a lattice ver- 
sion of the spinless p + ip model, with a potential V{x) — 
ae~^ /2ct ^ turning the hard-wall ending at x = effectively 
into a smooth end. The parameter a tunes the length scale over 
which the potential is turned on. The case a corresponds 
to a hard-wall boundary. The prefactor a has the dimension 
of energy and determines the height of the potential. For the 
calculations shown in the figure, we chose a = 5/i. The re- 
sults of the figure show a clear exponential dependence on a 
for states on the wire end terminated by the smooth function 
V, allowing for energies of the subgap states that are signif- 
icantly below the (already small) estimates for a rectangular 
geometry with hard-wall boundary conditions. 

The origin of the anomalously small energy splittings 
shown in Fig.|7]is the smoothness of all terms in the Hamilto- 
nian H. If the p+ip wire is coupled to a normal-metal stub, as 
in Sec. [Ill] and the superconducting order parameter A jumps 
at the interface a; = 0, no reduction of the end-state energies is 
found, even if the normal-metal stub is terminated by a smooth 
potential. (This scenario is well described by the calculation 
of Sec. [Ill]) On the other hand, one finds a suppression very 
similar to that shown in Fig. [7] if the superconducting order 
parameter goes to zero smoothly at the interface. We refere to 
Ref i52 for a discussion of the effect of a smooth confinement 
in the semiconductor model. 



VI. CONCLUSION 

In this article, we have investigated fermionic subgap states 
localized near the end of a spinless p-wave superconducting 
wires for two terminations of the wire — a normal-metal stub 
and a smooth confining potential — and in the presence of 
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weak disorder. The three scenarios give qualitatively different 
estimates for the energies of the subgap states. However, they 
share the common feature that a wire with N transverse chan- 
nels with a width W that is much smaller than the supercon- 
ducting coherence length ^ has int {N/2) fermionic end states, 
all with energy much below the bulk excitation gap A. These 
states appear for the topological phase (which has a Majorana 
fermion at the wire's end), as well as for the non-topological 
phase (which does not). 

The appearance of low-energy fermionic end states poses 
an obstacle to the identification of Majorana fermions through 
a measurement of the tunneling density of states at the wire's 
end, unless the energy resolution of the experiment is good 
enough to resolve the splitting between the fermionic end 
states. The corresponding energy scale 

^min scales propor- 
tional to A/fcp'^ ^ ^^//^ in the most favorable scenario we 
considered (wire's end coupled to a small normal metal stub), 
which is the same dependence as the subgap states in a vor- 
tex core.^^ The important difference with the subgap states in 
a vortex core is, however, that the number of fermionic end 
states is limited, so that there exists a maximum energy e^ax, 
whereas no such maximum energy exists in a vortex. Other 
terminations, such as a rectangular end with or without dis- 
order, or a smooth confinement potential, give significantly 
smaller values for Emin, and, hence, lead to stricter require- 
ments for the energy resolution required to separate an even- 
tual Majorana state from fermionic end states. 

The recent experiments that reported the possible ob- 
servation of a Majorana fermion involved semiconductor 
nanowires with proximity-induced superconductivity.^^ Ef- 
fectively, the induced superconductivity in these wires is of 
spinless p-wave type. However, it should be emphasized that 
this does not imply that the effective description of such a 
semiconductor wire with N transverse channels is a p + ip 
model with the same number of transverse channels. Instead, 
only those channels in the semiconductor that are effectively 
spinless {i.e., spin polarized or helical, depending on the rel- 
ative strength of the applied magnetic field and the spin-orbit 
coupling) appear in the effective description in terms of a 
p + ip model. (This latter distinction was overlooked in Ref. 
[30l ) Typically, this number is smaller than the number of 
transverse channels in the semiconductor In particular, the 
nanowires of the experiments of Refs.l22landl23iare believed 
to be thin enough that they map to a single-channel p + ip 
model. Hence, we do not expect that the mechanism for the 
generation of fermionic end states we consider applies to those 
experiments. However, it will apply to nanowires with a larger 
diameter, which we thus expect to exhibit a clustering of low- 
energy fermionic states in the topologically trivial as well as 
the topologically nontrivial phases. In this context, it is impor- 
tant to note that the condition that <C C does not a priori 
prevent the applicability of our analysis to thicker wires, be- 
cause the effective pairing potential A may decrease with W 
for proximity-induced superconductivity in the limit of thick 
wires (see Ref. 17 for an example in which A cx W^^). 

We gratefully acknowledge discussions with Felix von Op- 
pen and Inane Adagedeli. This work is supported by the 
Alexander von Humboldt Foundation in the framework of the 



Alexander von Humboldt Professorship, endowed by the Fed- 
eral Ministry of Education and Research. 



Appendix A: Relationship between the p + ip and proximity 
coupled semi-conductor models 

A practical realization of a the p + ip model can be found 
in semiconducting nanowires with strong spin-orbit coupling, 
laterally coupled to a standard s-wave superconductor and 
subject to a Zeeman field. In the following we discuss the 
precise relationship between the models. A related discussion 
can also be found in Ref.lT4l 

In two dimensions, and without coupling to the supercon- 
ductor, the Hamiltonian for this system reads 



P 

2m 



fi + Ba^ + aayp^ - a'a^Py, (Al) 



where a and a' set the strength of the spin-orbit coupling and 
B > is the Zeeman energy of the applied magnetic field. 
In the limit of a narrow wire (width W much smaller than 
the coherence length ^ of the induced superconductivity), sub- 
gap states as well as the above-gap quasiparticle states have a 
vanishing expectation value of the transverse momentum ky, 
which allows us to treat the transverse spin obit term as a per- 
turbation, initially setting a' = 0. Without the term propor- 
tional to a' different transverse channels do not couple to each 
other and the eigenfunctions of the Hamiltoruan Hj<s are of the 
form 



' \ J Vw 

where the angle 6k is defined as 



1 e'^^sin^, (A2) 



W ' 



sm «fe 



ak 



B 



cos Uk 



and the corresponding energies are 



(A3) 



-k,n 



2m 2mW^ 



± a/-B2 + a2fc2. (A4) 



Upon laterally coupling the semiconductor wire to an s- 
wave superconductor, the excitations are described by the 
Bogoliubov-de Gennes Hamiltonian 



BdG 



Aay -H;, 



P 

2m 



a CTxPy 



Ba, 
Aa. 



aCTyPi 



y 'X, 



(A5) 



where t^, Ty, and are Pauli matrices in electron-hole 
space. Without the transverse spin-orbit coupling a', the 
Bogoliubov-de Gennes Hamiltonian has a chiral symmetry, 
TyHTy — —H. In the basis of the normal-state eigenfunctions 
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the Bogoliubov-de Gennes Hamiltonian (A5i takes the 

form 



HsdG — 



- l^OyT^ COS 6'fe + Acr^Tj. sin 

- a'py{a^ cos 6'fe + ay sin 6'fe). 



(A6) 



In the Hmit, when both A and the spin orbit energy are 
smaller than the Zeeman splitting, the s-wave pamng term 
proportional to ay is ineffective, and each transverse chan- 
nel separately maps to two spinless p-wave superconductors, 
one for ip^ and one for V',7fc- Neglecting \ak\ in compar- 
ison to B, the corresponding pairing term A sin 6'fe(T2Ta; sa 
—A'kazTx, with 



A' 



aA 
~B~' 



(A7) 



Without the term proportional to a', the transverse channels 



in Eq. ( A6 1 can be treated independently (at least in the bulk of 
the wire, see below). If /i < S, only the "— " channels (eigen- 
spinors of az with eigenvalue —1) in Eq. ( |A6[ ) are topologi- 
cally nontrivial and can possibly have end states.^ Projecting 
the Bogoliubov-de Gennes Hamiltonian in the rotated basis 
( |A6| l onto these channels, one finds an effective Hamiltonian 
of the form 

iJ^ff = + - n- B]tz 

+ A'hkr^ + a'py , (A8) 
Without the transverse spin-orbit coupling a', the effective 



states at each end of the wire. The chiral symmetry is bro- 
ken by the transverse spin-orbit coupling a' . Because of the 
particle-hole symmetry of the Majorana modes, the matrix el- 
ements of this perturbation between the N Majorana end-state 
of with a' — are the same as the matrix elements of 

the p-wave superconducting pairing Hy of Eq. if we iden- 
tify A' — a' in the expression for Hy. 

If the condition /i < S is not met, the relation between 
the semiconductor and p + ip models is more complicated. 
For transverse channels for which h^ir'^n? /2mW'^ < fi — B 
the wire ends represent a chiral-symmetry-preserving pertur- 
bation that gaps out eventual Majorana end states, so that such 
channels may be disregarded when considering low-energy 
end states. For transverse channels for which 



(A9) 



Hamiltonian ( A8 1 has chiral symmetry and N Majorana end 



the Majorana end state in the "— band" (eigenspinors of az 
at eigenvalue —1 in the rotated basis) is protected in the pres- 
ence of the chiral symmetry, and only perturbations that lift 
the chiral symmetry can lead to a splitting of these end states. 
Projecting the Bogoliubov-de Gennes Hamiltonian in the ro- 
tated basis ( |A6| l onto these channels, one again an effective 
Hamiltonian of the form ( |A8| l, but with the additional restric- 
tion that only those transverse channels that meet the condi- 
tion ( |A9[ ) are considered. The number N of transverse chan- 
nels that meet this condition may be smaller than the original 
number of propagating channels in the semiconductor. 
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